Two-dimensional charge order in layered 2-1-4 perovskite oxides 
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Monte Carlo simulations are performed on the three-dimensional (3D) Ising model with the 2-1- 
4 layered perovskite structure as a minimal model for checkerboard charge ordering phenomena in 
layered perovskite oxides. Due to the interlayer frustration, only 2D long-range order emerges with 
a finite correlation length along the c axis. Critical exponents of the transition change continuously 
as a function of the interlayer coupling constant. The interlayer long-range Coulomb interaction 
decays exponentially and is negligible even between the second-neighbor layers. Instead, monoclinic 
distortion of a tetragonal unit cell lifts the macroscopic degeneracy to induce a 3D charge ordering. 
The dimensionality of the charge order in Lao.5Sri.sMn04 is discussed from this viewpoint. 

PACS numbers: 71.45.Lr, 05.50.+q, 64.60.Fr, 75.40.Mg 



Charge ordering (CO) phenomenon often appears in 
transition metal oxides including manganites, nickelates 
and probably cuprates. Of particular interest is the CO 
in layered perovskite oxides with K2MF4 lattice struc- 
ture (so-called 2-1-4 structure) such as Laa— x Sr x M O4 
(M=Ni,Mn,Cu) [1-5]. When the carrier concentration is 
a rational number, a commensurate CO can occur. Par- 
ticularly, in the case ofx = 1/2 or 3/2, the CO within the 
ab plane emerges with the checkerboard pattern at the 
wave vector Q — (n/a, tt /a), where a is a lattice constant 
within the plane. 

Neutron scattering measurements have revealed that 
below the transition temperature Tco = 217 K, 
Lao.sSr! 5M11O4 shows the checkerboard-type CO as a 2D 
long-range order (LRO) with a finite correlation length 
along the c axis [1] . This provides an evidence of a Bragg 
rod. In contrast, recent X-ray [6] and neutron [7] scat- 
tering experiments have indicated a Bragg peak at a 3D 
wave vector Q = (tt/o, n/a, 0) at a much lower tempera- 
ture than Tco, indicating a 3D CO [7]. The dimension- 
ality of this CO in the material remains controversial. 

To resolve this controversy, the frustration of the inter- 
layer Coulomb interaction shown in Fig. 1 plays a central 
role: The two relative configurations of the neighboring 
2D CO are completely degenerate when only the nearest 
neighbor interlayer coupling is taken into account. Frus- 
trated interaction yields (i) complicated patterns of the 
ordering such as an incommensurate state, (ii) suppres- 
sion of a LRO to realize the (quantum) liquid, (iii) glassy 
state, and so on [8,9]. The degeneracy is usually lifted in 
a nontrivial way. Even a simple order emerges due to the 
so-called 'order by disorder' mechanism [10], where the 
entropy of the fluctuation around each degenerate con- 
figuration differs and the system picks up the state with 



the largest entropy. 

Particularly for layered systems, when interlayer inter- 
action suffers from a frustration, the dimensional reduc- 
tion of a LRO occasionally occurs with a macroscopic 
degeneracy and only a weak universality relation is sat- 
isfied [11,8]: In special 3D models stacked with vertex 
models, a 2D LRO appears and critical exponents con- 
tinuously vary with an interlayer coupling constant [11], 
as in the eight- vertex model [12,13]. Therefore it is highly 
nontrivial and important to study what happens to the 
CO in the 2-1-4 structure in the presence of the interlayer 
frustration, which we undertake in this paper. 

We model the CO in terms of the Ising model by as- 
signing up and down spin configurations to Mn 3+ and 
Mn 4+ [14]. We consider the nearest-neighbor antiferro- 
magnetic (AF) interaction in the 2-1-4 lattice structure 
shown in Fig. 1 (a). The Hamiltonian is described as 

Here a\ = ±1 is the Ising variable defined at the site i in 
the £th layer. We have introduced 2D displacement vec- 
tors for nearest-neighbor sites r\ = (a, 0) or (0, a) within 
the plane and 8 = (a/2, a/2) between the adjacent layers. 
Hereafter, we consider the case of J > J± . In the limit of 
the independent 2D AF Ising models (J± = 0), a phase 
transition takes place at T 2D = 2/(sinh~ 1 1) ~ 2.269 [15]. 
Since J± is completely frustrated between the 2D AF or- 
dered planes, there occurs a 2 L --io\d degeneracy corre- 
sponding to two-fold for each layer (Fig. 1 (b)) where L z 
is the number of layers. 
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To study thermodynamic properties of the model (1) 
in detail, we perform Monte Carlo (MC) calculations. 
We employ the histogram algorithm [16] to obtain high- 
precision data, in addition to the metropolis algorithm 
with local flip as well as global or Wolff-type cluster flip 
[16]. We take a = 1 as a length unit and J = 1 as an 
energy unit hereafter. 

We expect that the AF order in each plane survives 
even with J± > 0. This is supported by MC results ex- 
emplified in Fig. 2 in the case of J± = 0.5: A temperature 
dependence of the specific heat exhibits a singularity at 
T ~ 2.18. At almost the same temperature, there occurs 
a systematic crossing of the Binder parameters g [16] for 
the 2D antiferromagnctism at the different system sizes, 
namely, g = 1 - (M 4 )/3(M 2 ) 2 with the 2D AF order 
parameter for a certain layer £, M = ^2iCr^c 1< ^' t . The 
bracket denotes the thermal average for the canonical en- 
semble. These consistently indicate that at least the 2D 
AF LRO within each plane emerges at a critical temper- 
ature T c ~ 2.18 which is slightly reduced from T 2B . 

Usually, a 3D LRO due to the order-by-disorder mech- 
anism [10] is likely to occur because the entropy force 
drives a 3D ordering in the presence of a finite inter- 
layer coupling. Despite this expectation, we found that 
there occurs no 3D ordering at all. Figure 3 shows 
the 3D spin structure factor S(k) at temperatures be- 
low T c estimated from Fig. 2. Here 5(k) is defined as 

S( k ) = £«' Ei^f^f V^""^ /^ 3 where k and 
Yu = (i, t) are 3D vectors. Since the weight is con- 
centrated on S(ir, 7r, k z ), we show the k z dependence of 
S(ir,ir,k z )/L 2 . There is no significant structure corre- 
sponding to the development of 3D Bragg peak. The 
inset shows a system-size dependence of S(ir,ir,k z )/L 3 , 
which clearly indicates that S(w, n, k z ) oc L 2 , i.e., the 2D 
LRO exists at every k z while there is no 3D ordering in 
the thermodynamic limit at all. Thus, there is a phase 
transition where the 2D AF correlation length diverges 
although the c-axis correlation length remains finite. 

This dimensional reduction is also confirmed by the 
low-temperature expansion of the free energy density 
with respect to J± [17], / = /o + f 2 + f± + ■ ■ ■ with 
/„ - 0((J±x.2D(Q)) n )- h and h are calculated 

as f 2 = 2X7^ fc (4J-iCOS^cos^) 2 X2D(fc) 2 and f* = 



expansion respect to J± is justified. Therefore, there is 
no possibility of the 3D ordering at any wave vector. 

Next, to understand the universality of this phase 
transition, we examine the critical exponents v and 
i] in the case of J± = 0.5. For this purpose, we 
perform MC simulations based on the histogram algo- 
rithm [16], to obtain high-precision data on the fol- 
lowing quantities in the vicinity of T c ; the 2D spin 
structure factor S 2 u(Q) = (M 2 )/L 2 , the 2D spin 
susceptibility X2d(Q) = {(M - (M)) 2 ), (d/d/3)\ogg, 
(d/d/3) log(M) and {d/d/3) log(M 2 }. We apply the finite- 
size scaling analysis to all the above quantities with 
S 2 d(Q)L^ 2 = ^(sL 1 ^), X2D (Q)L"- 2 - MsL 1 ^), 
[(d/d/3) log g\L z ~ x l v = ^(sL 1 '"), [(d/d/3) log(M) = 
$ 4 (ei 1/l/ ) and [(d/d/3) log(M 2 } = ^(eL 1 /"). Here, $'s 
are universal functions and e = (T — T c )/T c is a reduced 
temperature. We determine T c , v, rj and $i's from the 
least-square fittings of all the data at different L and T 
to $i with i = 1, • • • , 5. Figure 4 shows that this finite- 
size-scaling analysis works quite well without any loga- 
rithmic correction. The analysis gives T c = 2.177953(7), 
v = 0.8687(1) and r] = 0.24575(6). Statistical errors in 
the last digit are shown in parentheses. This high preci- 
sion is comparable to the best one obtained for the con- 
ventional 3D fsing model [16]. The set of these values 
for critical exponents is not consistent with either the 2D 
Ising (t] = 1/4 and v = 1) [18] or the 3D Ising universality 
class (rj = 0.032 and v = 0.63) [19]. This suggests that 
the present phase transition does not belong to either the 
2D or the 3D Ising universality class. 

For various values of Ji, we perform a finite-size- 
scaling analysis of S 2 n(Q) based on the MC results of 
the metropolis algorithm [20]. Figure 5 summarizes the 
estimates of T c and the critical exponents as a function of 
J±: Increasing J± from 0, T c decreases continuously from 
T 2D and the exponent v decreases continuously from the 
2D Ising value v = 1, while the exponent rj remains sim- 
ilar to the 2D Ising value rj ~ 1/4. Such continuous 
change of critical exponents appears in the eight-vertex 
model [12,13] and Ashkin- Teller model [21] and the so- 
called sliding phase of a 3D XY model [22]. We also 
found that in the bilayer system, the continuously vary- 
ing v but with almost fixed rj is obtained (Fig. 5). The 
universality class of this 2D phase transition is not deter- 
Sfci ••• fe 4 (rii=i 4 4J± cos ^f- cos % i )x2D(fei, k 2 , fe 3 , fe 4 ) 2 niined by a fixed point but by a fixed line, which is also 
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respectively. Here, \2D is the susceptibility for the 2D 
AF Ising model and \ 2 u is a four-point susceptibility 
defined by x 2D = T ' A : V.,........ ; (n?=i e" ife ^) x 

(HUiVji) - £/< CT i M i) CT W ^3)^(4))] • £' re P" 
resents the summation over three different combinations 
of (1, 2, 3, 4). The geometrical factor cos(fca;/2) cos(fc y /2) 
appears at all orders in the expansion, which is a charac- 
teristic feature of the complete frustration of the present 
model. Furthermore, since \id and \ 2 d behave as 
exp(— A/T) where A is the gap of the order of J, the 



found in some 2D frustrated systems [11]. Namely, the 
interlaycr coupling between the adjacent layers modifies 
the universality class of the model. 

To understand the dimensionality of the checkerboard 
CO experimentally observed in Lao.5Sri. 5 Mn04, we dis- 
cuss here the stability of the 2D LRO against possible 
3D LRO's due to (I) the long-range Coulomb interaction 
between layers and (II) lattice distortion. 

(I) First, we add to the model (1) the following inter- 
layer coupling term between the second-neighbors (Fig. 1 
(c)) to mimic the longer-range Coulomb interaction; 
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H' = with dP=Y.i° 



bility to J' ± is given by \4 = (C^ }/£ 2 - If X4 is finite 
for L — > oo, the 2D order is stable at least for small as 
in low-dimensional ordered phases in liquid crystals [23] . 
On the contrary our MC results for \i show a divergence 
as L 2 in the limit of L — > oo. This means that the spins 
within each layer act as a big spin of the size of L 2 . Then, 
the present 2D LRO is unstable against a J^-driven 3D 
LRO. However, this effect of J' ± is almost cancelled out 
by the remaining part of 1/r longer-range Coulomb in- 
teraction than J' ± between different layers: We define the 

average spin in the £th layer as = L~ 2 J2i a i^ e 1 ® 1 
and an effective interaction J^ ff (I — £') between and 
S^'\ We note jf(2£ + 1) = 0. There remains the 
complete frustration between the odd-spaced layers. Us- 
ing the Ewald summation, J^{21) is found to decay as 
e -4.5£c/a f or c j a Then even for the second-neighbor 
layers, J^ ff (2)/J is reduced to 10 
stants a = 3.86A and c = 12.442 
small energy scale of the order of /iK should be irrelevant 
in realistic experimental situations. 

(II) In realistic materials, there exists the electron- 
lattice coupling. In the present model in Eq. (1), if 
the lattice is deformed, the couplings J and/or Jj_ arc 
modified. There arc many candidates for the lattice 
distortion to lift the macroscopic degeneracy accompa- 
nied with the 2D LRO. Up to the linear order in dis- 
tortion, however, the system is unstable against only a 
monoclinic distortion of the tetragonal unit cell (Fig. 1 
(d)). This distortion produces two different interlayer 
couplings so that the frustration is partly removed. 
It is useful to introduce a structure factor correspond- 
ing to this distortion as F = (O^/L 3 , where On = 



The suscepti- 



8 with the lattice con- 
[24]. This extremely 
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with 5^ = (1/2, ±1/2) is the operator that directly cou- 
ples to the distortion. In the case of J± = 0.5, from 
the finite-size scaling of MC results shown in Fig. 6, we 
found that F diverges toward T c in the thermodynamic 
limit. The 2D order is unstable against the monoclinic 
distortion to lead a 3D LRO. If we rotate the unit cell by 
45° around the c axis, this distortion represents an or- 
thorhombic one. This is achieved by applying a uniaxial 
pressure along the a' or b' axis in the rotated frame. 

In conclusion, we have studied the phase transition and 
the dimensionality of the LRO in the AF Ising model 
on the 2-1-4 lattice: If the crystal structure has an ideal 
tetragonal symmetry, the 2D LRO appears with the non- 
trivial critical phenomena at T c . On the other hand, 
the system is unstable towards the monoclinic distortion, 
which replaces the 2D LRO with a 3D LRO. This con- 
clusion applies to the checkerboard CO phenomenon in 
Lao.5Sri. 5 Mn04. In the former case without the distor- 
tion, one should observe a stronger singularity in the spe- 
cific heat with the exponent a = 2(1 — v) >0 than the 



logarithmic divergence. In the distorted case, a change of 
the lattice constants b Q ^ 2a Q in the orthorhombic unit 
cell should be observed as an additional feature to a CO- 
induccd doubling of a unit cell and possibly a proposed 
pattern of distortion related to orbital order [25]. 
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FIG. 3. (Color) The 3D spin structure factor S(tt, 71", k z ) di- 
vided by the square dimension L 2 is plotted against k z below 
T c in the case of L = 64. The inset shows that S(w, w, 0)/L 3 
linearly goes to zero with 1/L. 




FIG. 1. (Color) (a) K2MF4 tetragonal lattice structure and 
the interactions in the model (1). The interlayer coupling 
J± is partially drawn for the clarity of the figure, (b) Frus- 
trated stacking of two adjacent planes. The solid and gray 
lattices are for neighboring planes projected to the stacking 
direction. Open and filled circles denote up and down spins, 
respectively, (c) Second-neighbor interlayer couplings (blue 
solid line) and the longer-range interactions between the sec- 
ond-neighbor layers (light-blue dashed lines), (d) Monoclinic 
lattice distortion to which the system shows an instability. 
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FIG. 2. Temperature dependences of (a) the specific heat 
and (b) the Binder parameter in the case of J± = 0.5. The 
lines are guides to the eye. 
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FIG. 4. (Color) The best-fit result of the finite-size scal- 
ing for correlation functions obtained by our MC calculation 
based on the histogram algorithm in the case of Jj_ = 0.5. 
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FIG. 5. (Color) Jj_ dependences of (a) the critical temper- 
ature, (b) the critical exponents v and r\. Dashed (dotted) 
lines denote the values of the 2D (3D) Ising model. Filled 
(open) symbols are for L 3 (L 2 x 2, namely bilayer) systems 
in model (1). Errorbars are omitted when they are smaller 
than the symbol sizes. Note that the values of Jj_ for 2-layer 
systems are divided by factor of 2 for convenience. 
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FIG. 6. (Color) The finite-size scaling for F obtained by 
our MC calculation in the case of J± = 0.5. We have used 
the values of T c and v obtained from the analysis given in 
Fig. 4. 
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